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Abstract 



We apply the realistic shell model which includes the coupling between many-particle 
(quasi-)bound states and the continuum of one-particle scattering states, called the 
Shell Model Embedded in the Continuum, to the spectroscopy of mirror nuclei: ^''F 
and ^"^O, as well as to the description of low energy cross section (the astrophysical S 
factor) for El, E2 and Ml components in the capture reaction: ^^0(p, 7)^''F. With 
the same microscopic input we calculate the phase shifts and differential cross-section 
for the elastic scattering of protons on ^^O. 
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I. INTRODUCTION 



Microscopic description of weakly bound exotic nuclei close to the drip-lines such as, 
e.g., or ^^Li in their ground state (g.s.), or nuclei close to the /9-stability line in the 
excited configurations like, e.g., ^^F in the first excited state J'" = 1/2^*", is an exciting 
theoretical challenge due to the proximity of particle continuum. Closeness of the scattering 
continuum implies that virtual excitations to continuum states cannot be neglected as they 
modify the effective interactions and cause the large spatial dimension of nuclear density 
distribution and, in particular, the existence of halo structures. Even basic concepts of the 
nuclear collective model hardly apply for those weakly bound configurations as the particle 
motion in the weakly bound orbit is presumably decoupled from the rest of the system [|l|. 
The microscopic description of such exotic configurations should treat properly the residual 
coupling of discrete configurations and the scattering continuum. Recently, it was proposed 
to approach this difficult problem in the quantum open system formalism which does not 
separate the subspaces of (quasi-) bound (the Q subspace) and scattering (the P subspace) 
states HJ^. Such a formalism can be provided by the Continuum Shell Model (CSM) 0] 
which in the restricted space of configurations generated using the finite-depth potential has 
been studied for the giant resonances and for the radiative capture reactions probing the 
microscopic structure of these resonances 

In this context, capture reaction ^^0(p, 7)^'^F is interesting for several reasons. First of 
all, precise experimental data in the energy range from 200 keV to 3750 keV are now available 
and the decays to the g.s. (J'^ = 5/2+) and to the first excited state {J"" = 1/2+) of ^^F 
have been accurately resolved 0. The strikingly different behavior of the astrophysical 
^-factor for the proton capture into the 5/2+ state and into the 1/2+ weakly bound state 
{Q = 105 keV) has been explained by the existence of proton halo in the 1/2+ state of ^^F 
0. Different theoretical approaches, including the potential model 0,|[, the model based on 
the Generator Coordinate Method (GCM) 0, or the K- and i?-matrix analysis [l^, have 
been tried to describe this reaction. Moreover, an expected simplicity of the low energy 
wave functions in ^^O, ^"^O and ^''F, allows to test certain salient features of models such as 
the effective interactions between Q and P subspaces and the possible quenching of matrix 
elements of the residual coupling which depends on the model space used in the calculations. 

Secondly, the exact knowledge of rate for the reaction ^^0(p, 7)^^F is necessary for mod- 
elling of nucleosynthesis process in the hydrogen-burning stars. Explosive hydrogen burning 
is believed to occur at various sites in the Universe, including novas. X-ray bursts, or the 



supermassive stars |]TT],|12|. Hydrogen burning of second-generation stars proceeds mainly 
through proton - proton chain and CNO cycle. The changeover from the pp chain to the 
CNO cycle happens near T ~ 2 ■ lO^K. The reaction ^^0{p, 7)^^F is of particular interest in 
this context as it provides a link to the higher branches of the CNO cycle. In particular, it 
starts the sidebranch : ^^0{p,'jY'^F{p,'jy^Ne{P+uy^F{p, aY^O. The contribution of CNO 
cycles to the total amount of produced energy in the Sun is small and CNO neutrinos ac- 
count only for about 0.02 of the total neutrino flux |TB[. Moreover, most of them are coming 



from the decay of ^^N and ^^O in the main CNO cycle CNO-I. Flux of '^^F-neutrinos', which 
is again two orders of magnitude smaller than the flux of neutrinos from ^^N and ^^O in 
CNO-I reactions, is controlled by the reaction ^^0(p, 7)^^F in the CNO-II cycle. Hence, the 
measurement of CNO neutrinos coming from different sources provides an accurate handle 
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on the thermonuclear reaction process in stars hke the Sun and allows in principle to distinct 
between different branches of the CNO cycle. 

The Shell Model Embedded in the Continuum (SMEC) model, in which realistic Shell 
Model (SM) solutions for (quasi-)bound states are coupled to the one-particle scattering 
continuum, is a recent development of the Continuum Shell Model (CSM) [^-0 for the de- 
scription of complicated low energy excitations of weakly bound nuclei. The SMEC approach 
is based on the realistic SM which is used to generate the A^-particle wave functions. This 
deliberate choice implies that the couphng between SM states and the one-particle scatter- 
ing continuum has to be given by the residual nucleon - nucleon interaction. In SMEC, 
like in the CSM, the bound (interior) states together with its environment of asymptotic 
scattering channels form a quantum closed system. Using the projection operator technique, 
one separates the P subspace of asymptotic channels from the Q subspace of many-body 
localized states which are build up by the bound single-particle (s.p.) wave functions and by 
the s.p. resonance wave functions. P subspace is assumed to contain {N — l)-particle states 
with nucleons on bound s.p. orbits and one nucleon in the scattering state. Also the s.p. 
resonance wave functions outside of the cutoff radius Rcut are included in the P subspace. 
The resonance wave functions for r < Rcut are included in the Q subspace. The wave 
functions in Q and P are then properly renormalized in order to ensure the orthogonality 
of wave functions in both subspaces. The application of the SMEC model for the descrip- 
tion of structure for mirror nuclei: '^B, ^Li, and capture cross sections for mirror reactions: 
^Be(p, 7)^B, ^Li(n, 7)^Li has been published recently 

We aim in SMEC at a microscopic description of low lying, complicated many-body 
bound and resonance states. For that reason, the description of particle continuum is re- 
stricted to the subset of one-nucleon decay channels. This should be a reasonable approach 
for describing the structure of mirror nuclei ^^F and ^''O around ^^O, and the capture reac- 
tions: ^^0(p, 7)^''F and ^^0(n, 7)^'^0. One expects that at higher excitation energies, e.g., 
above a emission threshold, the one-particle continuum approximation is too restrictive and 
the residual correlations generated in bound state wave functions by the coupling to those 
channels cannot be described. Effects of such correlations have been seen in the structure 
of J'^ = 3"*" resonances in ^B and ®Li which strongly couples to the three-particle decay 
channels. In this case, one should try to employ methods based on the cluster expansion 
of the wave function and the three-body continuum models M-jTBI . More complicated two- 



particle channels like, e.g., the a - decay channel, can be treated in SMEC, following the 
approach of Balashov et al. [1T7 . 



The paper is organized as follows. In Sect. II, we remind certain elements of the SMEC 
formalism and, in particular, those features of the S'-matrix in SMEC which are involved 
in the calculation of the elastic cross-sections, phase shifts and the capture cross-sections. 
Sect. Ill is devoted to the discussion of specific properties of ^^F and ^^O. In particular, 
properties of the matrix elements of the effective operator which couple Q and P subspaces in 
^^F and ^''O will be discussed in Sect. III.B. Features of the self-consistent average potentials 
in Q subspace for different many-body states of ^^F will be presented in Sect. III.C. Sect. IV 
is devoted to the discussion of ^^0{p, ■yY'^F capture cross sections for different multipolarities 
and different final states of ^''F. The differential elastic cross-section and the elastic phase- 
shifts in p -|- ^^O scattering will be compared to the experimental data in Sect. V. Summary 
and outlook will be given in Sect. VI. 
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II. THE FORMALISM OF SHELL MODEL EMBEDDED IN THE CONTINUUM 



The full solution of SMEC approach is constructed in three steps. In the first step, one 
calculates the (quasi-) bound many-body states in Q subspace. For that one solves the 
multiconfigurational SM problem : 



(1) 



where Hqq = QHQ is the SM effective Hamiltonian which is appropriate for the SM 
configuration space used. For solving (0), we use the code ANTOINE [0 which employs 
the Lanczos algorithm and allows for the diagonalization in large model spaces. 

For the coupling between bound and scattering states around ^^O, we use either a com- 
bination of Wigner and Bartlett forces (WB force) : 



Vi2 = -VS\a + (1 - «)Pr#(ri - ra) 



or the density dependent interaction (DDSMl) 



12 



P?2 in 



p{r)vZ{l - p{rW,lTyT, [p{r)v^- + (1 - p{r))vli)\ 5{r, - r^) 
is the spin exchange operator and p(r) in 

r — Tq' 



IS : 



p(r) 



1 -|- exp 



d 



(2) 



(3) 



(4) 



with ro = 2.64 fm and d = 0.58 fm. The DDSMl force (§) has four parameters determined 

115.5, all in units MeV-fm^. Similar 



as: V, 



00 



36.45, 



^00 



-297.9, 



109.4 and 



force has been used by Schwesinger and Wambach [|T^ for the description of giant resonances 
in (Ip - Ih) - {2p - 2h) space. However, as compared to their interaction which is close to 
the Landau - Migdal type of interactions, the DDSMl parameters f). 



00' "^00' '^oi 



Vm and Dgf , are 



reduced by a constant factor 0.68. Also the parameter ro has been somewhat reduced to 
better fit the experimental matter radius in oxygen. For the WB residual interaction (^, 
the overall strength Vi2^ is adjusted to reproduce as good as possible the energy spectra and 
the decay widths of states in nuclei around ^^O. A reasonable compromise is provided by 
^12'* = 300 MeV-fm'^. The relative contribution of spin exchange term has been discussed 
for the p-shell nuclei *B and ^Li [Q. It was argued that the spin exchange part is very 
small ((1 — a) ~ 0) for most p-shell nuclei but is expected to increase for heavier nuclei. 
In present studies, we use (1 — a) = 0.27, which is a standard value in the IsOd mass 
region [Q. It should be stressed also that matrix elements of Hqq are not modified by the 
residual coupling between P and Q because they are fitted to experimental discrete levels 
and narrow resonances in the continuum and, therefore, they are believed to contain already 
those coupling effects. 

The SM wave function has an incorrect asymptotic radial behavior for unbound states. 
Therefore, to generate both the radial s.p. wave functions in the Q subspace and the scatter- 
ing wave functions in P subspace we use the average potential of SW type with the spin-orbit 
part included: 
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U(r) = Vofir) + Vso^li2l-s)-^ + Vc , (5) 

where = 2fm^ is the pion Compton wavelength and /(r) is the spherically symmetrical 
SW formfactor : 

(6) 

The Coulomb potential Vc in (H) is calculated for the uniformly charged sphere with radius 

Rq. 

For the continuum part, one solves the coupled channel equations: 

{E^^^ - Hpp)it^ ^ Y.iE^^^ - H^c')&^ = , (7) 

c 

where index c denotes different channels and Hpp = PHP . The superscript (+) means 
that boundary conditions for incoming wave in the channel c and outgoing scattering waves 
in all channels are used. The channel states are defined by coupling of one nucleon in the 
scattering continuum to the many-body SM state in (A^ — l)-nucleus. The channel - channel 
coupling potential in (|^) is : 

H^, = {T + U)5,,+vi. , (8) 

where T is the kinetic-energy operator and vj^^ is the channel-channel coupling generated by 
the residual interaction. The potential for channel c in (^) consists of 'initial guess', U{r), 
and diagonal part of coupling potential vj.^ which depends on both the s.p. orbit (pij and the 
considered many-body state J'^ . This modification of the initial potential U{r) change the 
generated s.p. wave functions (pij defining Q subspace which in turn modify the diagonal 
part of the residual force, etc. This means that the solution of coupled channel equations 
(0) is accompanied by the self-consistent iterative procedure which yields for each channel 
independently the new self- consistent potential : 

[/(-^)(r) = U{r)+vi^'''\r) , (9) 

and consistent with it the new renormalized matrix elements of the coupling force. The 
parameters of the initial potential U{r) are chosen in such a way that U'^^^^r) reproduces 
energies of experimental s.p. states, whenever their identification is possible. 

The third system of equations in SMEC consists of inhomogeneous coupled channel 
equations: 

- Hpp)u^^ = HpQ<^, = w, (10) 

with the source term Wi which is primarily given by the structure of - particle SM wave 
function $j . The explicit form of this source is given in [Q. These equations define functions 
djl^"* , which describe the decay of quasi-bound state $j in the continuum. The source 
Wi couples the wave function of A^-nucleon localized states with (A^ — l)-nucleon localized 
states plus one nucleon in the continuum. Formfactor of the source term is given by the 
self-consistently determined s.p. wave functions. 



fir) 



exp 



r-Kn 
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The full solution of SMEC equations is expressed by three functions: , and uji 

m ■■ 



where 



H'S = Hqq + HqpG^^^Hpq , (12) 



is the effective SM Hamiltonian which includes the coupling to the continuum, and G^p^ is 
the Green function for the motion of s.p. in the P subspace. Matrix Hqq is non-Hermitian 
(the complex, symmetric matrix) for energies above the particle emission threshold and 
Hermitian (real) for lower energies. The eigenvalues Ei — are complex for decaying 
states and depend on the energy E of particle in the continuum. The energy and width 
of resonance states are determined by the condition: Ei{E) = E [|]. The eigenstates 
corresponding to these eigenvalues can be obtained by the orthogonal but in general non- 
unitary transformation : 

^•. = EA,<f, (13) 

j 

where the coefficients Pij form a complex matrix of eigenvectors in the SM basis satisfying : 

f^ikPjk = Sij (14) 

k 



Analogously, one can define : 



and 



EA.c^, (15) 



= l>i + cDi . (16) 

Inserting them in (|TT|), one obtains : 

n = + E^^T^^r^TTT^ < I Hqp I Ce > (17) 
i E- Ei + {i/2)Ti 

for the continuum many-body wave function projected on channel c, where 

= + E dE ^(^) _ ^, < I HpQ I > , (18) 

is the wave function of discrete state modified by the coupling to the continuum states. It 
should be stressed that the SMEC formalism is fully symmetric in treating the continuum 
and bound state parts of the solution : represents the continuum state modified by 
the discrete states, and fij represents the discrete state modified by the coupling to the 
continuum. 
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A. Radiative capture and electromagnetic transitions 



Having obtained the full solution of SMEC, one can calculate many observable quantities. 
E.g., the calculation of capture cross-section goes as follows. The initial SMEC wave function 
for the [p(8>(^ ~ 1)]"'' ' system is : 



vi/,(r) = y: 



X xT X X 



It 



i-h) 



(19) 



and the final SMEC wave function for the A^-system in the many-body state JV is : 



bSJb 



kjb 



{Jf) 

nif 



(20) 



It and s denote the spin of target nucleus and incoming proton, respectively, ^ifsj^^ is the 



coefficient of fractional parentage and Ui^j^ is the s.p. wave in the many-particle state Jf . 
These SMEC wave functions (0,^) are then used to calculate the transition amplitudes 
T^^ and T*^^ for El , E2 and Ml transitions, respectively 0,^. The radiative capture 
cross section is : 



a 



~9~\K — — ■ 



he J \hcj 2s + 1 2It + 



E1,M1 1 2 



(21) 



a 



R2 



75 \klJ \hcj \hcj 2s + 1 2/t + 1 ^ ' 



E2 1 2 



(22) 



and fi stands for the reduced mass of the system. 

Another interesting quantities are B{EX), B{MX) transition matrix elements between 
the SMEC wave functions J^' and Jj^ for the [A^] - system : 



i<vi/K^;oiiQ(nA)i|vi/.(jr)> 

2J, + 1 



(23) 



where 11 = E (electric) or M (magnetic), and A = 1,2... (multipolarity). O(nA) in 
the electromagnetic transition operator . 



IS 



B. Properties of the S'-matrix 



The S'-matrix for the scattering of a nucleon by a nucleus is given by the asymptotic 
behavior of the total wave function (iTi 



E-E, + {il2)V 



(24) 
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with the incoming wave only in the channel Cq. The asymptotic conditions for the solution 
(P^) have been analyzed in Ref. 0]. Here we only give the final result for the amplitude of 
the partial width : 

71") = - exp(-.5(°)) (^^) Y.PmY.j dr^f(r)wf(r) , (25) 

where 5^^^ is the background scattering phase : 

5f = ^arg^L°^ , (26) 

and S^^^ is determined from the asymptotic, large distance behavior of Using the 

proportionality relation between matrix elements in ( p^ and the amplitudes of partial width 
(p5|), one derives the S"- matrix elements : 

^cco = S^l - ^ exp[.(5(°) + 5^)] ^ . (27) 

Amplitudes of the partial widths 7^") as well as the real En and imaginary r„ parts of 
complex eigenvalues of the effective shell model Hamiltonian Hqq (|T^) entering the S- 
matrix elements are explicitly energy dependent. 7^"'' are also complex due to the channel - 
channel coupling and the mixing of quasi-bound states embedded in the continuum . The 
partial width for channel c can then be defined by : 

rn,c=|7^f , (28) 



though the total width r„ is no longer sum of partial widths 



El7l"f = rnE 1/5-1' >f. (29) 

c i 

because the eigenvectors are normalized in a sense of (0), what implies : 

El/5™|'>1 . (30) 



III. SPECTROSCOPY OF ^^F AND ^^O NUCLEI 

A. The effective interactions for the Oplsld SM space 

For the purpose of the present study we have constructed SM effective interactions in 
the cross-shell model space connecting the Op and IsOd shells. The interactions have three 



distinctive parts: the Op-shell part taken to be the CK(8-16) interaction |22], the IsOd-shell 
part taken to be the Brown- Wildenthal interaction ||23|. For the cross-shell matrix elements, 
we use the G matrix of Kahana, Lee and Scott (KLS) [|24[. This part of the interaction has 



to be modified phenomenologically. We proceed as follow: the IsOci s.p. energies above the 
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^He core are taken from generalized monopoles fit of energies of single-particle and single- 
hole states over the mass table 



25| . Then the cross monopoles are adjusted to reproduce 



^^N, ^^O and ^'^O, ^^F spectra. The center of mass components are removed with the usual 
Gloeckner - Lawson prescription [^] adding a center of mass kinetic energy hamiltonian. For 



natural parity states, we consider only Ohu excitations and for non natural parity states, Ifno 
excitations. Of course, it is well known and has been studied in the past, that ^^O (and^^O, 
^^F) have significant amount of 2p — 2h {3p — 2h) and 4p — Ah {5p — 4h) admixtures in their 
low lying states Nevertheless, the shell model calculation of these admixtures is not 

straightforward: 



• Firstly, it has been already stressed in the past by different authors p9|-pT[], that Nhuj 
mixing converges slowly with N and, if stopped for example after N = 2, causes a 
strong lowering of Op — Oh states due to 2hu} pairing correlations which were already 
present in pure 2p — 2h states. One could overpass this artifact of the calculation by 
artificially lowering the Nhu unperturbed configurations through the monopoles of 
the interaction, but the control of convergence of the mixing would be impossible as 
higher order admixtures would destroy the picture. 

• Secondly, the building blocks of our interaction are made of OHlj Op and IsOd phe- 
nomenological interactions, which in principle already contain, in particular through 
pairing renormalizations, the higher order correlations. For this reason, we have de- 
cided to fix the unperturbed (4p — 4/i) 0+ state around its experimental energy (6.049 
MeV). 

Again, the derivation of precise correlated wave functions for ^^O and ^''O, ^^F is out of 
the scope of this paper and we will refer to the mixing derived by other studies P^,P 



B. The effective operator of the residual coupHng 

Expressions for the matrix elements of residual interaction which couples P and Q sub- 
spaces have been given in Refs. @,@]. There are two kinds of coupling operators. The channel 
- channel coupling reduced matrix elements involve one-body operators : 

O^s = i<^y^sf (31) 

Diagonal part of O induce the renormalization of the s.p. average potential. Matrix elements 
of the operators O are calculated between different many-body states in the (A^ — l)-system, 
i.e., in ^^O, and they depend sensibly on the amount of 2p — 2h and 4p — 4h correlations in 
the g.s. of 1^0. 

Reduced matrix elements of the source term Wi in the inhomogeneous coupled channel 
equations (|I0|) contain a product of two annihilation operators and one creation operator : 

KsiDp = {4{~a,~as)y- . (32) 

Operators TZ enter in the calculation of complex eigenvalues of Hqq. Matrix elements of 
TZ are calculated between different initial state wave functions in either [p{2)(A^ — 1)] or 
[n0(A^ — 1)] system, i.e., in ^''F or ^^O, and a given final state wave function in ^^O. These 
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matrix elements depend on the configuration mixing in the many-body wave functions of 
i^F, and 1^0. 

We perform the SM calculations for ^^O ground state in which admixture of 2p — 2h and 
4p — Ah configurations are neglected : 

|i6o, 0+ > = aooKOp - Ohy='> > + a22\{2p - 2hy='> > + a^H^p - 4/i)^=° > + ■ ■ ■ 

~aoo|(Op-0/i)^=° > (33) 

i.e., the g.s. of ^^O is taken to be a pure {Op^^) Op — Oh configuration and the second excited 
state with predominantly 4p — Ah structure is not mixed with the ground state. How big 
are these np — nh admixture depends on the SM effective interaction and the effective s.p. 
space used p7| , p8| . In Table H we give results of different SM calculations for 0^^ in ^^O. 



Likewise in ^^O and ^^F, the g.s. J'^ = 5/2^ and the excited states = 1/2^, 3/2^ 
correspond in our approximation to the pure configuration of one particle in, respectively, 
0(i5/2, lsi/2 or 0(^3/2 shells outside of the (Op^^)-core of ^^O, i.e.: 

|^^F;J^7^ = + > =ai;^\{lp-Ohy'' > + ai;\{3p-2hy^ > + ■■■ 

~af;|(lp-0/i)-'' > . (34) 



The amount of these admixture calculated by Brown-Green |23] and Zuker-Buck-McGrory 
p8| can be read from the Table |1[ The experimental spectroscopic amplitudes for 5/2i and 
1/25*" states can be found in Table |T|. 

The negative parity states involve the excitations across the shell closure and involve one 
hole in the p shell and two particles in the IsOd shells. In our model of ^^O where 2p — 2h 
and 4p — Ah configurations are neglected, the spectroscopic factors of negative parity SM 
states are equal zero. Hence the particle decay width for these states is coming from the 
functions Ui (|T^) which describe the continuation of Q space wave functions into P, i.e., 
from the modification of discrete states by the coupling to the continuum (18). 

Neglecting higher order correlations in the g.s. of ^^O and in positive parity excited states 
of either ^^F or ^^O means that residual coupling between Q and P which involves those 
states should be properly rescaled, i.e., the matrix elements of O, TZ should be quenched : 

< ^^O, J/|C|^^0, Ji >il^^^^hed)^ ^^^^Y ^ 16q^ J/|C|^^0, Ji >i^on-guenched) ^35^ 

< 1^0, Jf\n\'^0, Ji >(9"^"^'^^'^)= aooflio < '^O, J/|7^|^^0, Ji >i"on-quenched) ^gg^ 

< ^^O, Jf\n\'^F, Ji >(l^-^-hed)^ ^^^^A ^ 16q^ J/|7^|^^F, Ji >i^on-quenched) ^37^ 

Amount of configuration mixing for negative parity states is large (two particles in IsOd 
shell and 1 hole in Op shell) and therefore matrix elements of TZ are not quenched. In our 
calculations, we take the values given by Brown - Green ||2^ for the amplitudes aoo and af^ 
(see Tables I and ID . 
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C. The self-consistent average potentials 



Construction of the Q subspace in SMEC is achieved by the self-consistent, iterative 
procedure which for a given initial average s.p. potential (^) and for a given residual two- 
body interaction between Q and P yields the self-consistent s.p. potential U^^^\r) which 
depends on the s.p. wave function ipij, the total spin J of the A^-nucleon system as well as 
on the one-body matrix elements of (A^ — 1) - nucleon daughter system. As explained above, 
the positive parity states in ^''F and ^''O are described as a pure one particle configuration 
in IsOfi shells. Consequently, the spectroscopic factors for the states J'^ = 1/2^ ,5/2'*' 
and 3/2"^ are equal 1. Consistently with this approximation, we can identify position of 
proton (neutron) s.p. orbits lsi/2, 0^5/2 and 0^3/2 with, respectively, J'^ = 1/2'*', 5/2^ 
and 3/2'*' many body states of ^''F (""^"^O), i.e., we may ask that U^^'^\J^) provides the 
energy of s.p. orbit at the energy of the corresponding many-body state with respect to the 
threshold. This choice is essential for quantitative description of radiative capture cross- 
section. Consequently, the initial SW potentials U{.P) have all different depth parameters. 
Moreover, they depend on the choice of the residual interaction coupling P and Q. On the 
contrary, we keep a common radius Rq = 3.214 fm, a common diffuseness a = 0.58 fm and 
a common spin-orbit strength Vso = 3.683 MeV for all of them. The depth parameters 
of initial SW potentials for [pl^^^O] and [n{g)^^0] systems, are summarized in Table [rV| . 
For a given nucleus (^''F or ^^O), all those different central potentials correspond to the 
one and the same equivalent potential f/^^^^(r) in the SW form which binds s.p. orbits at 
the same energy as obtained in the self-consistent potentials for neutrons and protons. For 
both WB and DDSMl residual interactions, the depth of f/(^«) is Vq = -52.46 MeV in ^^F 
and Vq = —52.49 MeV in ^''O. The potential radius, the diffuseness and the spin-orbit 
strength are the same as used for the initial potentials U{J'^). The spectrum of single- 
particle (-hole) energies in this potential is given in Table |V| for ^^F and in Table |V]] for 
^^O. These energies are determined from the data of ^''F, ^''O (for particles) and ^^N, ^^O 
(for holes). The difference of s.p. energies ep(lsi/2) — ep(0c?5/2) = 0.495 MeV in ^"^F and 
en(l'Si/2) — en(0ci5/2) = 0.871 MeV, which should be compared with the excitation energy 
of = l/2i state in those nuclei, is a manifestation of the Thomas-Ehrman shift |3S 



and is included effectively in our equivalent and self- consistent s.p. potentials for all studied 
channels. 

In Fig. |l]we show examples of calculated potentials in ^^F for the proton s.p. orbitals lsi/2 
and 0(i5/2 in the total spin states: J'^ = 1/2+ and J'^ = 5/2+, respectively. The calculations 
have been performed using the appropriate initial potentials U{.P) (see Table |V[) for the 
DDSMl residual interaction (|^). For example, potential U (5/2+) is chosen in such a way that 
the self-consistent potential f/(''^)(5/2+) (see Table [V|) yields Oo?5/2 proton s.p. orbit bound 
at the experimental binding energy of the ground state (g.s.) = 5/2^. Similarly, the 
choice of f/(l/2+) and the determination of U'^^'^\l/2'^) is associated with the reproduction 
of the experimental binding energy of the first excited state J'^ = 1/2^. The J'^-independent 
equivalent potential yields the 0c?5/2 and lsi/2 proton s.p. orbits at the position given for 
them in the corresponding U^^^\5/2~^) and U^^^\l/2'^) potentials. 

The self-consistent potential (the solid line) strongly deviates from the SW form. In 
the center U^^'^^r) has a strong maximum which is absent in the initial potential U{r) 
(the dashed line). The self-consistent potentials [/("^)(l/2+) and [/("'^)(5/2+) are different 
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(compare the solid curves on l.h.s. and r.h.s. of Fig. p, in spite of the fact that up to spin- 
orbit couphng the equivalent potential f/(^'J) (the dotted lines) in these states is identical. 
This clearly shows how strong is the state dependence of both the self-consistent average 
fields and the renormalized matrix elements of the coupling force. One should also notice the 
difference in surface region between t/("^)(l/2+), U^'^\h/2+) and f/^"'?). In (5/2+) one 
may notice decrease of the potential radius in comparison with the radius of U'^^'^\ However, 
as compared to the initial average potentials ^7(5/2"'") and [/(1/2"'"), both the self-consistent 
potentials [/^*'^-*(5/2+), U'^^'^\l/2'^)^ and the equivalent potential f/^*^"?) are deeper and their 
effective radii are bigger. 

In Fig. we show examples of calculated potentials in ^^F for the proton s.p. orbital 
Opi/2 in the many-body states J'^ = 1/2^. The calculations have been performed using the 
initial potentials U{l/2-) for the WB (l.h.s. of the plot) and DDSMl (r.h.s. of the plot) 
residual interactions. For both interactions, the equivalent potential U'^^'^\r) is the same. 
It is interesting to notice how different are the self-consistent [/(*'^) (r) potentials for the two 
considered residual interactions. As a rule, the renormalization of initial s.p. potential is 
weaker for the density-dependent DDSMl interaction. Similarly as for the self-consistent 
potentials in J'^ = 5/2'*' and J'" = 1/2^ states (see Fig. ), the effective radius of U^^^\r) 
shrinks as compared to the radius of U^^'^\r) . Again this effect is stronger for the WB 
residual interaction. 

In general, the surface region of average potential shows weak sensitivity to the self- 
consistent correction. In our case, the radial dependence of self-consistent correction for 
all positive parity states in ^'^F and ^''O is given by the Opi/2 and 0p3/2 radial formfactors 
corresponding to the well bound single-hole state and not by the weakly bound single-particle 
states lsi/2 or 0^5/2. For that reason, even for the = 1/2^ state in ^"^F which is bound 
by 105 keV, induced renormalization of the surface in the self-consistent potential decreases 
the radius of potential in the surface region. 

D. Discussion 

In this section, we shall present the SMEC results for the spectrum of mirror nuclei: 
^^F and ^^O. The results depend mainly on: (i) the nucleon - nucleon interaction in Q 
subspace, (ii) the residual coupling of Q and P subspaces, (iii) the self-consistent average 
s.p. potential which generates the radial formfactor for s.p. bound wave functions and s.p. 
resonances, and (iv) the cutoff radius for s.p. resonances. Below, we shall comment on their 
relative importance. 

1. Spectrum of^'^F 

Fig. 1^ compares the SM energy spectrum of ^^F for positive parity (l.h.s. of the plot) 
and negative parity (r.h.s. of the plot) states, with those obtained in the SMEC approach 
for WB and DDSMl residual interactions. The experimental data are plotted separately for 
positive and negative parity states as well. SM calculation is performed using the effective 
interaction described in Sect. III. A. In the column denoted 'SMEC (WB)' we show results of 
SMEC approach with the residual coupling between Q and P subspaces which is given by the 
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mixture of Wigner and Bartlett forces with the spin-exchange parameter (1 — a) = 0.27. 
The overall strength parameter Vi2^ = 300 MeV-fm^ has been obtained by fitting the spectra 
of and ^^O. In the column denoted 'SMEC (DDSMl)', the DDSMl density dependent 
residual interaction (|]) is used. Only the coupling matrix elements between the J'^ = Of g.s. 
wave function of ^^O and all considered states in ^^F are included. Zero of the energy scale 
for ^^F is fixed at the experimental position of J'^ = 1/2]^ first excited state with respect 
to the proton emission threshold for all different examples of SMEC calculations which are 
shown in Figs. ^ and 0. 

The iterative procedure to correct U{r) and to include the diagonal part contribution 
of residual interaction has been described in the previous section. The self-consistently 
determined s.p. potential U^^'^\r) is then used to calculate radial formfactors of coupling 
matrix elements and s.p. wave functions. Different initial potentials U{J'^) (see Table [TV] ) 
are used for the calculation of self-consistent average potentials for different many-body 
states in ^^F. These different self-consistent potentials correspond to a unique equivalent 
SW potential U^^'^\r). The renormalization of initial potential by the residual coupling of Q 
and P subspaces is the same for all states J'^ of the same spin and parity. For positive parity 
states, as a result of the restriction in SM calculations in Q subspace, the renormalization 
of average potential operates only for s.p. orbits which have the same spin j and parity vr as 
those of the many body states J, vr. E.g., for the many body state J'^ = 5/2"*" of ^'^F (^"^O), 
only potential of the 0(^5/2 s.p. orbit is modified by the coupling to the continuum. Similarly 
for = 1/2+ many body state and the s.p. orbit lsi/2. For those s.p. orbits which in a given 
many body state are not modified by the selfconsistent renormalization, we calculate radial 
formfactors using the 'universal' equivalent s.p. potential U^^'^\r) which does not depend 
neither on the many body state nor on the s.p. orbit. Of course, this special property is only 
due to a simple structure of ^^O and ^^F nuclei in our SM calculations. For neutrons in ^^F 
there is no renormalization of the average potential and we use the equivalent potential for 
neutrons to get the radial formfactors (see Table 0). Supplementary informations concerning 
the results shown in Fig. |^ can be found in Table |V11| . 

The spectrum of ^''F is insensitive to certain approximations in the SMEC. G.s. energy 
relative to the proton emission threshold is reasonably well reproduced by the SMEC. Cou- 
pling to the continuum induces strong relative shifts of 5/25*" and 3/2^ states with respect to 
the 1/2^ state and the negative parity states. The position of g.s. with respect to the energy 
threshold for proton emission changes by about 500 keV due to the inclusion of coupling 
to the g.s. of ^^O (compare columns denoted 'SM' with those denoted 'SMEC (WB)' and 
'SMEC (DDSMl)' in Fig. ^ and in Table |V11| ). These shifts are generally smaller for the 
DDSMl residual coupling. The coupling to excited 0^ state at E* = 6.049 MeV in ^^O (0^ 
has predominantly Ap — 4h structure with respect to the dominant Op — Oh configuration in 
the g.s. 0^ ) is expected to be unimportant in ^^F or ^^O. Calculated width of states depends 
on chosen residual couplings (compare 'SMEC (WB)' with 'SMEC (DDSMl)') and in gen- 
eral the agreement is better for the density dependent interaction. The effective quenching 
of the coupling matrix elements (^5 H37|) for positive parity states does not solve this problem 
completely. To get a better agreement with the data, the energy separation of 5/2^ and 
1/2^ SM states should be decreased while the splitting of 5/25*" and 3/2^ SM states should 
remain unchanged. This means that the separation of Is subshell and the centroid od Od 
subshells should be increased by few hundred keV. We did not however pursue the studies 
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in this direction and we keep standard values of s.p. energies p5| 



The effect of quenching of the couphng operator between Q and P can be seen in Fig. ^ 
where the comparison between the SMEC calculations with the quenching factors (entry 
'SMEC (DDSMl)') and those without quenching factors (entry 'SMEC (DDSMl*') are 
shown. The strong downshift due to the continuum coupling is seen for the (5/2 j*" - 3/2]'") - 
pair with respect to both the negative parity states and the 1/2^ state which moves down- 
wards much less. The effect of quenching of the coupling operator on negative parity states 
is, as expected, very weak. 

A useful measure of the radial wave function is the B{E2) transition matrix element 
between l/2f and 5/2f bound states. In SMEC, we obtain 70.3 e^fm^ for DDSMl residual 
coupling. This value has been calculated using the effective charges : Cp = 1.41 and e„ = 0.47 
and the BG quenching factors (Table The experimental value for this transition is 
B{E2)e^p = 64.92 e^fml The SM prediction for this transition is B{E2)sm = 32.7 e^fm^ 
assuming the oscillator length b = A^^^ fm and the same effective charges as used in the 
SMEC calculations. The better agreement of calculated in SMEC B{E2) values with the 
data is mainly due to more realistic radial dependence of the lsi/2 s.p. orbit in J'^ = l/2f 
many body state. The rms radius for this orbit in SMEC is < >^^^ = 5.212 fm, as 
compared to the value of < >^^^ = 3.629 fm for the 0^5/2 s.p. orbit in = 5/2+ state. 
Nuclear rms radius of ^^F is 2.664 fm in the g.s. and 2.814 fm in the first excited state 1/2"*". 



2. Spectrum of O 

Fig. 1^ compares the SM energy spectrum of ^''O for positive parity (l.h.s. of the plot) 
and negative parity (r.h.s. of the plot) states, with those obtained in the SMEC approach 
for WB and DDSMl residual interactions. The experimental data are plotted separately 
for positive and negative parity states as well. SM calculation is performed using the same 
effective interaction as used for ^''F spectrum shown in Figs. |^ and In the column denoted 
'SMEC (WB)' we show results of SMEC approach with the residual coupling given by the 
mixture of Wigner and Bartlett forces (0). In the column denoted 'SMEC (DDSMl)', the 
density dependent interaction DDSMl is used as the residual interaction. Parameters of the 
residual couplings are the same as used to describe the ^^F spectrum. Only the coupling 
matrix elements between the J'" = Of g.s. wave function of ^^O and all considered states in 
^^O are included. For all different examples of SMEC calculations which are shown in Fig. |^, 
zero on the energy scale is fixed at the experimental position of J'^ = 5/2^" state which is 
close to the neutron emission threshold. 

The self-consistently determined s.p. potential U^^'^\r) is then used to calculate radial 
formfactors of coupling matrix elements and s.p. wave functions. Different initial potentials 
U{J'") for [n (g) ^^O] (see Table [I^) are used for the calculation of self-consistent average 
potentials for different many-body states in ^^O. These different self-consistent potentials 
correspond to a unique equivalent SW potential U'^^'^\r) (see Table [V^ . For protons there 
is no correction from the residual interaction modifying the average potential and we use 
the equivalent potential for protons to get the radial formfactors (see Table [V|) . 

The energy intervals between the g.s. and the positive parity states l/2f and 3/2^ in 
SMEC(DDSMl) calculations perfectly reproduce the experimental energy sequence for these 
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levels. The WB residual force leads to a too strong relative shift of 5/2^ with respect to 
l/2f and 3/2f states. The quenching factors have been introduced in the calculations for 
positive parity states similarly as discussed before for ^^F. One should remind that in ^^F, in 
spite of using the quenching factors, the 5/2f state was shifted downwards with respect to 
other positive parity states. Improved situation in ^^O is due to the larger separation of this 
state with respect to the particle continuum. Nevertheless, the coupling to the continuum 
is necessary here to bring closer the 3/2f and l/2i SM states. The shift of the ground state 
energy in SMEC relative to the neutron emission threshold is solely due to the choice of zero 
on the energy scale at the experimental position of 5/2^ weakly bound state. Calculated 
width of states depends on chosen residual couplings (compare 'SMEC (WB)' with 'SMEC 
(DDSMl)') and in general the agreement is better for the density dependent interaction. 

As a measure of the radial wave function in ^^O we calculate the B{E2) transition 
matrix element between 1/21" and 5/2f bound states. In the SMEC, we obtain S.le^fm^ 
for DDSMl residual coupling. Changing the residual coupling and using WB interactions 
leaves this value almost unchanged. The B{E2) value in SMEC has been calculated using 
the effective charges: Cp = 1.41 and e„ = 0.47, and the BG quenching factors (Table y). The 
experimental value for this transition is B{E2)f.xp = 6.2e^fm^ whereas the SM prediction 
is B{E2)sM = 3.62 e^fm^ assuming the oscillator length h = A^/^im. and taking the same 
effective charges as in the SMEC calculations. The close agreement of calculated in SMEC 
B{E2) values with the data is due to a more realistic radial dependence of the lsi/2 s.p. 

orbit in J'" = l/2i many body state. The rms radius for this orbit is < >^^'^ = 4.184fm, 
as compared to the value: < >^^^ = 3.439 fm for the 0^5/2 s.p. orbit in Jf = 5/2"^ many 
body state. 

IV. THE ASTROPHYSICAL FACTOR FOR ^^0{P,-fy'^F 

In Fig. ^ we show the calculated multipole contributions {El, E2, Ml) to the total 
capture cross section as a function of the cm. energy Ecm, separately for the transitions to 
the g.s. 5/2^ and to the first excited state 1/2^ in ^^F. The SMEC calculation is done with 
the DDSMl residual interaction as used in the calculations of spectra shown in Figs. 3 and 
4. Parameters of initial potentials U{5/2~^) and U{l/2~^) can be read from Table [rV| . The 
proton threshold energy is adjusted to agree energies of calculated and experimental l/2i 
state with respect to the proton emission threshold. The photon energy is then given by 
the difference of cm. energy of [^^O + p] j. system and the experimental energy of the final 
state Jf in ^'^F. 

The dominant contribution to the total capture cross-section for both 5/21" and 1/2^ 
final states, is provided by El transitions from the incoming p wave to the bound 0^5/2 and 
lsi/2 states. We took into account all possible El, E2, and Ml transitions from incoming 
s, p, d, f, and g waves but only El from incoming p - waves give important contributions. 
In the transition to the g.s., the El contribution from incoming fY/2 wave gives is by a 
factor ~100 smaller than the contribution from P3/2 wave for energies of incoming proton up 
to 3.5 MeV. The E2 contribution is smaller by at least three orders of magnitude in both 
branches. 

In the transition to the first excited state 1/2^, the dominant Ml contribution comes 
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from the incoming si/2 wave and is by ~2 orders of magnitude smaller than the El contri- 
bution. The Ml contribution is important, in particular, in the branch ^^0(p, 7)^'^F( J'^ = 
5/2]^) for EcM > 3MeV, where its contribution is of a similar order as the El contribu- 
tion due to 3/2^ resonance in the ^3/2 scattering wave. One should also notice that the 
energy dependence of Ml component in the capture cross section is totally different in 
^^0(p,7)^^F(J^ = 1/2^) and ^^0(p, 7)^^F( = 5/2^) branches. The contribution of Ml 
for ^^0(p, 7)^^F( J'^ = l/2f) is decreasing with increasing energy above Eqm ~ 1 MeV and 
approaches zero for Eqm ~ 3.5 MeV. The Ml contribution begins to grow again at higher 
energies. Since a mi cross section for 1/2^ final state originates only from the absorption of 
s-wave proton and, moreover, since by construction: < \E'/|\I^.j > = 1 (the spectroscopic fac- 
tor of 1/2^ state equals 1), therefore the energy behavior of cta/i cross section is governed by 
the Breit-Wigner like factor coming from the 1/2^ many-body wave function (see eq. ( pUj ) 
or (p!7[) ). So (Jmi depends crucially on the energy dependence of this single eigenvalue. While 
the real part of E{J'^ = l/2i) eigenvalue behaves smoothly with energy, the imaginary part, 
on the contrary, after initial rise becomes almost equal to zero at Ecm — 3.5 MeV and starts 
to rise again for higher energies. Since we have only one J'^ = 1/2+ state in our SM space, 
the energy behavior of its eigenvalue determines directly the behavior of a mi cross-section, 
because in the limit: F ^ 0, the Breit-Wigner factor behaves like the 5-function. In the 
more complete SM calculations which yield more J'" = 1/2"*" states, this unusual effect is 
expected to be reduced. Nevertheless, because its principal cause is the single-partial-wave 
characteristic of this transition, we believe that the trace of it remains. 

In the branch ^^0{p, 7)^^F( J'^ = 1/2^), one may also notice small resonant part in El at 
the position of 1/2 resonance, but otherwise the cross-section has non-resonant behavior. 
The non-resonant El contribution is a good measure of spatial extension of 5/2^ and 1/2^ 
wave functions respectively. For the decay branch ^^0(p, 7)^'^F( J'^ = 1/2^'"), due to the 
very small binding energy of J'^ = l/2i state {Q = 105 keV) and its particularly simple 
structure, this cross section is a sensible measure of the extension of the proton lsi/2 orbit. 
It is essential for the calculated cross-section that the lsi/2 proton orbit in the self-consistent 
average potential for J'^ = 1/2"*" is bound by 105 keV. Even small modification of this value 
by different choice of the depth parameter Vq in f/(r), introduce the modification of ^ei 
which is larger than any modification due to possible uncertainties in the potential radius 
Ro or its surface diffuseness a. 

Fig. shows the total S'-factor as a function of the cm. energy, as well as its values for 
the ^^0(p,7)^'^F(J'' = l/2i) and ^^0(p, 7)^^F( J'^ = 5/2|) branches. The overall agreement 
with the experimental data of Morlock et al is very good. The resonance contributions of 
l/2i in ^^0{p, 7)^'^F( J'^ = 1/2^^) is well reproduced even though their precise energy position 
is slightly shifted. This resonance contributes only through its P subspace continuation 
because the spectroscopic factor for this state is equal zero. The contribution of 5/2^ 
resonance in the ^^0(p, 7)^^F( J'^ = 5/2f) branch is so narrow in energy that has been 
omitted in the figure. Energy of this resonance in SMEC (DDSMl) is Eqm = 3.44 MeV. 

The strong increase of the S'-factor at low energies for 1/2+ component is a direct 
consequence of the proximity of continuum for lsi/2 s.p. orbit in = l/2i state. The 
energy dependence of S'-factor as Ecm — ^ has been fitted by a second order polyno- 
mial to theoretical points obtained in the interval from 13 to 50 keV in steps of 1 keV. 
For the DDSMl residual couphng we have S'(O) = 9.25 x 10"^ MeV-b, and the loga- 
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rithmic derivative is S'{0)/S{0) = —4.98 MeV ^. The ratio of Ml to El transition is 
5mi(0)/^ei(0) = 7.55 x lO^'l The El, Ml, E2 contributions from different incoming waves 
at EcM can be seen in Table The E2 and Ml components of the ^-factor for 
the branches ^^0{p,'jy^F{J'^ = 1/2^) and ^^0{p,'-fy^F{J'^ = 5/2f) summed over all partial 
incoming wave contributions are shown in Fig. For the WB residual coupling, El and E2 
contributions are only ~ 0.7% smaller but Ml is about 20 times bigger than for the DDSMl 
case, so the ratio of Ml to El transition is 5'mi(0)/S'ei(0) = 0.154. The ratio of E2 and 
El contributions for ^'^0(p, 7)i^F( = 1/2^) is: a^^/a^^ = 1.62 x 10~^ 2.22 x 10"^ and 
5.425 X 10"^ at 20, 100 and 500 keV, respectively. 

For the deexcitation to the g.s. 5/2^, the fit of calculated S'-factor as Ecm yields: 
S{0) = 2.89 X 10"^ MeV-b and S'{0) /S{0) = 0.767 MeV^^ The ratio of Ml to El transition 
is 5'mi(0)/5'ei(0) = 3.1 x 10~^. El, Ml, E2 contributions from different incoming waves at 
Ecm can be seen in Table 0. It is particularly interesting to notice large E2 contribution 
from si/2 incoming wave having large negative logarithmic derivative. The ratio of E2 and 
El contributions for i60(p, 7)i^F( J'^ = 5/2^) is: a^^/a^^ = 8.5 x lO^^, 3.5 x 10"^ and 
1.35 X 10~'^ at 20, 100 and 500 keV, respectively. In the total cross section this ratio is: 
a^'^/a^^ = 4.4 x 10-^ 3.7 x 10"^ and 6.1 x 10"^ at 20, 100 and 500 keV, respectively. 



V. ELASTIC CROSS-SECTION AND PHASE-SHIFTS 

Another observable quantities that can be calculated using the solutions of the SMEC are 
the elastic phase shifts and elastic cross-sections for different proton bombarding energies. 
Results for the elastic phase shifts, shown in Figs. and those for the elastic cross-section, 
shown in Fig. |T0|, have been obtained using the SMEC solutions for identical parameters 
of initial potentials and the DDSMl residual coupling between Q and P subspaces as dis- 
cussed before in Sect. III.D.l for the spectra of ^^F and in Sect. IV for the capture reaction 
^^0(p, 7)^''F( J''). The elastic phase shifts are well reproduced by the SMEC for all states 
except for the 5/2+. In this case the calculated and experimental |^| values differ by about 



5°. This discrepancy disappears when the resonant contribution to this phase shift is re- 
moved from the SMEC solution (compare the solid and dashed lines in the plot for 5/2"*" 
in Fig. It may be surprising at first sight to see so large resonant contribution to the 
phase shift associated with the g.s. wave function. For energies below the proton emission 
threshold the coupling to the continuum introduces only the hermitean modifications of the 
Hamiltonian which shift the energy of 5/2^ state in SMEC with respect to its initial position 
given by the SM but do not generate any width for this state. The coupling of Q and P 
subspaces is non-local and, hence, the effective Hamiltonian in SMEC is energy dependent. 
For excitation energies above the proton threshold, the coupling of Q and P subspaces gen- 
erates the non-hermitean corrections to the effective Hamiltonian which yield the imaginary 
part of the eigenvalue and generates the resonant like behavior which is so well visible in the 
5/2^ elastic phase shift. The origin of this large 'halo' of 5/2^ bound state in the continuum 
for positive energies and the large shift of the real part of the 5/2i eigenvalue for negative 
energies (see Fig. ^ and ^) is the same and should be traced back to the incorrect description 
of correlations in the SM wave function for this state. As discussed in the previous chapter, 
one expects that the 5/2^ state in ^''F and ^^O is not a pure s.p. excitation outside of the 
^^O core [27]. Therefore, it is also natural to quench the matrix elements oi Q - P coupling 
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as described in Sect. III.B. This procedure, which consists of using quenching factors for 
positive parity states to correct for missing higher order np — nh components in their wave 
function, is actually not sufficient as the example of elastic phase shift for 5/2f demonstrates 
clearly. The quenching procedure cures the problem of unphysically large energy shift for 
the 5/2i state due to the coupling to the continuum but does not correct the wave function 
for missing correlations which in turn lead to the disappearance of 'halo' of this state for 
positive energies. This deficiency of present SM calculations in ^^F and ^''O, which is seen 
consistently both in the energy spectrum as well as in the elastic phase shifts, demonstrates 
potentiality of SMEC approach for mass-regions far off the (3 - stability valley where the 
experimental information about exotic nuclei will be scarce and one will have to use both 
the spectroscopic information and the information from the scattering experiments to learn 
about the structure of those nuclei. This example shows also that in the SMEC approach 
which unifies description of discrete state properties and the scattering continuum, one may 
use different kinds of experimental data to fix those few parameters of the model such as 
the overall strength of the residual Q - P coupling or the radius and diffuseness of the initial 
average potential. 

Elastic excitation functions at a laboratory angle of 166° calculated in SMEC with 
DDSMl residual interaction are compared with the experimental data in Fig. |T^. The 
calculated values at low energies are too low as compared to the data and this is again due 
to the too strong halo of 5/2f state for positive energies (compare the dashed and solid 
lines in Fig. ^ at low energies). Removal of the resonant contribution from this state to 
the elastic cross section brings the calculations close to the data. The agreement between 
calculated and experimental low energy cross-sections provides a supplementary check of 
the spatial extension of self-consistent potential and, hence, of the radial form factors of s.p. 
wave functions. One should however keep in mind that, in general, the information from 
the elastic cross-sections may be strongly perturbed by the non-locality effects in the SMEC 
effective Hamiltonian which depend strongly on the many-body correlations in the SM wave 
functions as the above example demonstrates. 

On the average, the agreement between experimental and SMEC results for the elastic 
excitation functions is reasonable if one keeps in mind large sensitivity of this experimental 
measure to even small inaccuracies in the energy position and width of resonances. The 
model predicts correctly the interference pattern due to 1/2^ resonance at Ep ~ 2.6 MeV 
and 3/2^ resonance at Ep ~ 5.2 MeV. Also, the 5/2^ resonance is correctly predicted 
by the calculations though its width is too narrow to be presented in the figure. For the 
same reason, this resonance was not plotted in the proton capture cross-section for the 
"^^0(p, 7)^'^F( J'^ = 5/22*") branch in Fig. [7[ The difference between interference patterns in 
the data and in the SMEC calculation for 3.8 MeV < Ep < 5 MeV is mainly due to the 
reversed order of 3/2j~ and 3/2^ resonances in SMEC calculations as compared to the data. 
Consequently, the hole in the experimental elastic excitation function at Ep ~ 4.2 MeV, due 
to 3/2]^ resonance interfering with the 3/2^ resonance cuts sharply only the high energy tail 
of the 3/2^ resonance. 
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VI. SUMMARY AND OUTLOOK 



In this work we have apphed the SMEC approach for the microscopic description of ^^F 
and ^^O spectra, the low-energy radiative capture cross sections in the reaction ^^0{p, 7)^^F, 
and the elastic cross section for the reaction ^^0(p,p)^^0. In the SMEC model, which is 
a development of CSM model for the description of low energy properties of weakly 
bound nuclei, realistic SM solutions for (quasi-)bound states are coupled to the one-particle 
scattering continuum. For that reason, we use realistic SM effective interaction in the Q 
subspace and introduce residual force which couples Q and P subspaces. (For this residual 
coupling we take either a combination of Wigner and Bartlett forces or the density dependent 
DDSMl interaction which is similar to the Landau - Migdal type of interactions.) This 
deliberate choice of interactions implies that the finite-depth potential generating P space 
and matrix elements of the residual Q — P coupling, have to be determined self-consistently. 
The self-consistent iterative procedure yields then new state-dependent average potentials 
and consistent with them new renormalized matrix elements of the coupling force. These 
renormalized couplings and average potentials are then consistently used both in Q and P 
subspaces for the calculations of spectra, capture cross-sections, elastic cross-sections, elastic 
phase-shifts, etc. 

Simultaneous studies of spectroscopy in mirror systems as well as different reactions 
involving one nucleon in the continuum, allow for a better understanding of the role of dif- 
ferent approximations and parameters in the model. The dependence on radius, diffuseness 
or spin-orbit coupling parameters of the initial potential U{r) is not very important and 
they can be taken from any reasonable systematics. On the contrary, the depth of U{r) has 
to be carefully adjusted so that the energies of s.p. orbits in U^^^\r) for [n^{N — 1)] and 
[p0{N — 1)] systems, whenever their identification is possible, reflect the binding of many- 
body states near the particle emission threshold in the nucleus A^. This is very important for 
the quantitative description of reaction cross-sections. In the case of ^''F (^''O), correct iden- 
tification of lsi/2 and 0^5/2 s.p. orbits and hence the determination of an appropriate depth 
parameter in f/(l/2+) and f/(5/2+) is unambiguous because the spectroscopic amplitudes 
in l/2i and 5/2^ are close to 1. Different binding of mirror nuclei ^^F/^''0 leads for the 
same J'^ many-body states to different U'^'^^\r) in these nuclei. This in turn causes breaking 
of mirror symmetry of SM spectra for these nuclei and is, e.g., an essential ingredient in 
understanding the difference between B{E2) values for the transition 5/2f — > l/2f in ^^F 
and ^^O, as discussed in Sect. III. 

SMEC model in its present form includes the coupling to one-nucleon continuum. The 
wealth of experimental data can be described in a unified framework of SMEC in this 
approximation. These include: (i) the calculation of energy spectra, -B(nA) transition ma- 
trix elements and various static nuclear moments such as the magnetic or mass/charge 
quadrupole moments etc., (ii) the calculation of various radiative capture processes: (p, 7), 
(n, 7), Coulomb breakup processes: (7,p), (7, n) and elastic or inelastic cross sections {p,p), 
{n,n ); some of these observables have been discussed in this work. Problem of isospin sym- 
metry breaking due to the coupling to the continuum can be addressed by comparing elec- 
tromagnetic processes, e.g., B{IIX) transition matrix elements for certain states in mirror 
nuclei, and weak interaction processes like the first-forbidden /3-decay in mirror reactions. 
Finally, for nuclei close and beyond the proton (neutron) drip lines, the spontaneous proton 
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(neutron) radioactivity can be studied in the microscopic framework of SMEC (SM). These 
unifying features of SMEC approach are extremely useful for understanding of the structure 
of exotic nuclei far from the /3 - stability for which the available experimental information 
will be scarce. 

In this work we have studied nuclei close to the doubly magic ^^O in order to understand 
certain basic features of the SMEC and, in particular, of the Q - P coupling operator acting 
in the restricted SM configuration space. The resulting quenching of operators O and TZ ( 
Eqs. (|3T| ) and (|32|) respectively) could be related to the spectroscopic amplitudes for positive 
parity states in ^^F (^^O) and to the amount of 2p — 2h,4:p — Ah correlations in the g.s. of 
^^O. It was found also that this SM motivated correction of the effective operator does not 
solve the problem of 'halo' of discrete states for positive energies. This problem results from 
non-locality of the effective SMEC Hamiltonian and, more precisely, from the non-hermitean 
corrections to the eigenvalues for positive energies which generate the imaginary part. The 
imaginary part of eigenvalues and, hence, the size of this halo effect, is particularly large 
for pure single-particle (single-hole) configurations. For this reason, the simplification of 
structure of the many body states by neglecting the configuration mixing can in certain 
cases lead to an unphysical enhancement of resonant-like correction from bound states in, 
e.g., the elastic cross-section or the elastic phase-shifts. In this sense, ^''F and ^^O nuclei, 
with predominantly s.p. structure of positive parity states and extreme sensitivity to higher 
order correlations in the many-body wave functions are somewhat pathological. We believe 
that this problem will disappear in nuclei having more particles in the open shells, in which 
case SM will produce sufficient amount of mixing in the wave functions. 

More complicated decay channels involving, e.g., a particle, ^He or '^H in the continuum, 
are beyond the scope of SMEC in its present form. The future extension of the SMEC for 
such cluster configurations is possible in a framework proposed by Balashov et al |T^. It is 
encouraging, however, that these possible shortcomings in the description of decay channels, 
are so unambiguously reflected in the calculated decay width for these states . In general, 
the decay width is particularly sensitive to the details of the SM wave functions involved 
and to the values of matrix elements of residual coupling so they provide a sensible test of 
the quality of SMEC wave functions and/or approximations involved. 

The present studies have shown that SMEC results depend sensitively on very small 
number of parameters. Some of them, like the parametrization of the residual interaction 
which couples states in Q and P subspaces, has been established in the present work for 
S(i-shell nuclei. The others, related to the quenching of the effective coupling operator can 
be explained consistently with the SM analysis of spectroscopic amplitudes. Finally, the 
energies of s.p. states, which determine the radial wave function of many-body states, are 
bound by the SM spectroscopic factors and experimental binding energy in studied nuclei. 
This gives us a confidence that the SMEC can have large predictive power when applied to 
the nuclei in the less known regions of the mass table. The calculations can be performed on 
a similar level of sophistication as the SM which, with the recent progress in SM techniques 
and the effective interactions has been applied to the medium-heavy nuclei 
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TABLES 



TABLE L The comparison of amplitudes of SM wave functions in five different models for 
O. The entry WBT92 corresponds to the results of Warburton, Brown and Millener 



OJ^ state in 
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using the 'A4;i^ '-method |^. The entry WBP(4p4/i) corresponds to the results of SM calculations 
using WBP interaction |^ and the 'A4p4/i '-method ||3l|. Results of the Brown-Green (BG) p7t , 
Zuker-Buck-McGrory (ZBM) |^8| and Haxton-Johnson (HJ) are also given. 



Model 



BG 13] 
ZBM [H 



HJ H 
WBT92 15 
WBP(4p4/i) Ul\ 



{Op - Ohy^=^i 



0.874 
0.71 
0.648 
0.748 
0.775 



{2p - 2hy^=^t 



{Ap - Ahy^=^t 



0.469 
0.58 
0.67 
0.574 
0.557 



0.130 

0.14 
0.333 
0.299 



TABLE IL The amplitudes of dominant {Ip — Oh) component in the SM wave function for 
5/2 (the ground state) and excited states 1/2^, 3/2]*", calculated by Brown-Green (BG) ||2^ and 
Zuker-Buck-McGrory p^. 



Model 


{ip-ohy^='/^t 


{Ip - ohy^='^/^t 


{Ip - ohy^=^i^t 


BG [27| 
ZBM |28 


0.881 
0.65 


0.718 



0.901 
0.69 



TABLE in. Experimental spectroscopic amplitudes for the ground state ^/2^ and first excited 
state 1/2| 





Ref. i 


Ref. m 


Ref. il 


Ref. |3^ 


5/2+ 


0.949 


1.14 


0.964 


1 


1/2+ 


1 


0.866 


0.916 


0.905 
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TABLE IV. The parameters of initial potentials U{r) (^) used in the calculations of 
self-consistent potentials U^^'^\r) for the WB and DDSMl residual interactions. U'^^'^\r) are con- 
structed for various positive and negative parity states in ^^F and ^''O. For all considered cases 
the radius of the potential is Rq = 3.214 fm and the diffuseness parameter is a = 0.58 fm. The 
spin-orbit parameter is Vso = — 3.683 MeV. For more details, see the description in the text. 



bystem 


J [MeV J 


h 


Vq [iVieVJ (WB) 


Vq [MeVJ (DUhJvllj 




5/2^ 


f\ 1 

04/2 


—42.228 


—42.123 




1/2^ 


lSl/2 


—44.485 


—45.875 




o /o-l- 

3/2^ 


0^3/2 


—42.416 


—42.237 




1/2 


lPl/2 


—39.448 


—42.330 




3/2- 


lP3/2 


-39.914 


-42.200 




5/2- 
7/2- 


0/5/2 
0/7/2 






[n (g) 1^0] 


5/2+ 


0^5/2 


-42.060 


-42.123 




1/2+ 


l'Sl/2 


-44.588 


-45.875 




3/2+ 


0c^3/2 


-42.440 


-42.237 




1/2- 


lPl/2 


-41.310 


-42.330 




3/2- 


lP3/2 


-41.613 


-42.200 




5/2- 


0/5/2 


-43.921 


-42.213 




7/2- 


0/7/2 


-43.878 


-42.701 



TABLE V. The energies of s.p. orbits in the equivalent average potential U^^'^\r) for protons 
and neutrons in ^'^F which yield the same s.p. energies as the self-consistent potentials U^^'^\r; J'^). 



Ij 


Vq (Protons) 


Vq (Neutrons) 


OSl/2 


-28.168 




0P3/2 


-14.287 


-18.452 


Opi/2 


-12.131 


-16.297 


0(^5/2 


-0.6 




l'Sl/2 


-0.105 




0(^3/2 


2.497 




0/7/2 


12.27 




0/5/2 


22.287 
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TABLE VI. The energies of s.p. orbits in the equivalent average potential U^^'^\r) for neutrons 
and protons in ^"^O which yield the same s.p. energies as the self-consistent potentials U^^'^\r; J'^). 



Ij 


Vq (Neutrons) 


Vq (Protons) 




-32.823 




0p3/2 


-18.447 


-14.8 


Opi/2 


-16.378 


-12.73 


04/2 


-4.143 




l'Sl/2 


-3.273 




04/2 


-0.91 




0/7/2 


8.701 




0/5/2 


16.397 





TABLE VIL SM energies and SMEC energies and widths vs. experimental ones of ^''F nucleus. 
The proton separation energy is adjusted in order to reproduce the binding energy of the first 
excited 1/2"*" state. Different labels denote as follows: '(WB)' - results of SMEC calculations for 
the WB interaction (^) , '(DDSMl)' - results of SMEC calculations for the density dependent 
DDSMl interaction (^). In these two cases, the quenching factors (35 - 37) for the positive parity 
states have been included. Label (DDSMl*) denotes results of SMEC calculations without the 
quenching factors. Only the g.s. of ^^O is included in the coupling matrix elements. The cut-off 
radius is Rcut = 9.5 fm for the 4/2 s.p. wave function in 3/2^ many body state which is in the 
continuum. For the details of the residual interaction which couples Q and P subspaces, see the 
discussion in the text. 





EsM 
(MeV) 


E (WB) 
(MeV) 


r (WB) 
(keV) 


E (DDSMl) 
(MeV) 


r (DDSMl) 
(keV) 


E (DDSMl*) 
(MeV) 


r (DDSMl*) 
(keV) 


Eexp 

(MeV) 


r 

exp 

(keV) 


5 + 
2 


-0.826 


-1.467 




-1.377 




-1.675 




-0.600 




1 + 
2 


-0.105 


-0.105 




-0.105 




-0.105 




-0.105 




1- 
2 


2.134 


2.486 


15 


2.286 


4.4 


2.388 


6.2 


2.504 


19 


5- 
2 


3.279 


3.654 


0.04 


3.44 


~0 


3.544 


~0 


3.257 


1.5 


3 + 
2 


5.044 


4.085 


859 


4.096 


926 


2.876 


848 


4.400 


1530 


9- 
2 


3.786 


4.169 


~ 


3.946 


~ 


4.050 


~ 


4.620 




3- 
2 


4.459 


4.747 


152 


4.551 


120 


4.636 


154 


4.040 


225 


7- 
2 


5.045 


5.299 


9 


5.172 


1.2 


5.256 


2.2 


5.072 


40 


5- 
2 


5.016 


5.307 


1 


5.175 


0.02 


5.279 





5.082 


< 0.6 


1- 
2 


5.156 


5.538 


2.9 


5.313 


5 


5.416 


7 


5.437 


30 


3- 
2 


5.277 


5.661 


2.9 


5.432 


14 


5.535 


17 


4.888 


48 
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TABLE VIII. The same as in Table VII but for ^'O. For more details see the description in 
the text. 





EsM 


E (WB) 


r (WB) 


E (DDSMl) 


r (DDSMl) 


E 


r 

exp 




(MeV) 


(MeV) 


(keV) 


(MeV) 


(keV) 


(MeV) 


(keV) 


5 ' 
2 


-4.406 


-5.245 




-4.938 




-4.143 




-1 -4- 
2 


-3.685 


-3.855 




-3.856 




-3.273 




1 ~ 
2 


-1.446 


-1.464 




-1.450 




-1.088 




5 
2 


-0.301 


-0.301 




-0.301 




-0.301 




2 


1.464 


0.309 


83 


0.543 


208 


0.942 


96 


9- 
2 


0.206 


0.207 


~ 


0.206 


~ 


1.073 


< 0.1 


3- 
2 


0.879 


0.767 


49 


0.804 


40.4 


0.410 


40 


7- 
2 


1.465 


1.440 


0.08 


1.297 


0.4 


1.554 


3.4 


5- 
2 


1.436 


1.428 


~ 


1.396 


0.06 


1.589 


< 1 


1- 
2 


1.576 


1.575 


1.9 


1.572 


3.8 


1.796 


32 


3- 
2 


1.697 


1.699 


1.5 


1.691 


8.8 


1.236 
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TABLE IX. El, Ml, E2 contributions in the limit Ecm from different incoming waves to 
the 5-factor for the branch ^^0{p,j)^''^F{J'^ = 1/2^;''). The values are extracted by fitting second 
order polynomial to the calculated SMEC values for the DDSMl residual coupling. For more 
details see the description in the text. 



i60(p, 7)^'F(1") 


MC 


El 


Ml 


E2 




Pl/2 


P3/2 


Sl/2 


d3/2 


d5/2 


5(0) [MeV-b] 


3.06 X 10-3 


6.12 X 10-3 


6.93 X 10-^ 


5.40 X lO-'^ 


8.11 X lO-'^ 


S'{0)/S{0) [MeV-i] 


-4.95 


-4.95 


-9.45 


0.023 


0.025 



TABLE X. The same as Table ^ for the transitions in the branch ^^0{p,-fy^F{J'' = 5/2J) 







i60(p, 7)1^ 








MC 


El 


Ml 




P3/2 


hj2 


f7/2 


d3/2 


d5/2 


S{Q) [MeV-b] 


2.83 X 10-"^ 


1.11 X lO-'^ 


2.22 X 10-6 


1.43 X 10-*^ 


8.65 X lO-'^ 


S'(0)/S(0) [MeV-1] 


0.534 


7.85 


7.86 


0.36 


3.36 


MC 


E2 




Sl/2 


d3/2 


d5/2 


§7/2 


§9/2 


5(0) [MeV-b] 


3.46 X 10-6 


7.60 X 10-9 


3.26 X 10-*^ 


1.6 X 10-^2 


2.0 X 10-^^ 


S'(0)/S(0) [MeV-i] 


-17.76 


3.15 


3.28 


18.3 


18.3 
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FIGURES 




FIG. 1. The self-consistent potentials for lsi/2 s.p. orbit in J'^ = 1/2^ many body states (the 
l.h.s. picture) and 0^5/2 s.p. orbit in J'^ = 5/2"*^ many body states (the r.h.s. picture) in ^^F are 
calculated for the DDSMl residual interaction (^). 
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FIG. 2. The self-consistent potentials for Opi/2 s.p. orbit in J'^ = 1/2" many body states 
in ^"^F are calculated with the WB residual interaction (l.h.s. of the plot) and the DDSMl 
interaction (P) (r.h.s. of the plot). 
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FIG. 3. SM and SMEC in different approximations labelled 'SMEC (WB)', SMEC (DDSMl)' 
vs. experimental states of ^''F nucleus. The proton threshold energy is adjusted to reproduce 
position of the \/2'\^ first excited state. The shaded regions represent the width of resonance 
states. For the details of the calculation see the description in the text and in the caption of 
Table IVDI. 
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FIG. 4. SM and SMEC in different approximations labelled 'SMEC (DDSMl*)' (no quenching 
factors for positive parity states), SMEC (DDSMl)' (with quenching factors) vs. experimental 
states of ^''F nucleus. The proton threshold energy is adjusted to reproduce position of the l/2f 
first excited state. The shaded regions represent the width of resonance states. For the details of 
the calculation see the description in the text and in the caption of Table VII. 
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TT = + 



71 = - 



2 
1 

1 
2 
3 
4 



- neutron threshold 



Hiiiniiiiinin 3/2 



7/2 
3/2 
9/2 

3/2 
5/2 

1/2 



SM SMEC SMEC 

(WB) (DDSMl) 



1/2"" 
5/2^ 



EXP 



SM 



SMEC SMEC 
(WB) (DDSMl) 



EXP 



FIG. 5. The same as in Fig. |^ but for ^''O nucleus. The proton threshold energy is adjusted to 
reproduce position of the 5/2 j" excited state. The shaded regions represent the width of resonance 
states. For the details of the calculation see the description in the text and in the caption of 
Table IVlTl. 
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FIG. 6. Multipole contributions to the capture cross section to the g.s. ^^0(p, 7)^^F( J'^ = | ) 
and to the first excited state ^^0{p, 7)^''F( J'^ = ^~^) are plotted as a function of the center of mass 
energy Eqm- The SMEC calculations for El (the solid line), E2 (the dotted line) and Ml (the 
dashed line) have been done for the DDSMl residual interaction (^). For more details see the 
description in the text. 
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0.008 

^ 0.006 




^ 0.004 
co^ 0.002 



1 2 



FIG. 7. The astrophysical S'-factor for the reactions ^^0(p, 7)^'^F(J'^ = 5/2|,^) and 
^^0(p, 7)^''F( J'^ = 1/2^) is plotted as a function of the center of mass energy Ecu- The ex- 
perimental data are from [|^. The SMEC calculations (the solid line) have been done for the 
DDSMl residual interaction (H). For more details see the description in the text. 
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I I I I I I I I I 

0.5 1 1.5 2 

Er.. [MeV] 



FIG. 8. The E2 and Ml components of the astrophysical S-factor for the reactions 
i60(p,7)i^F(J'^ = 5/2+) (the dashed hnes) and ^''0(p, 7)i^F(J'^ = 1/2+) (the sohd Unes) is 
plotted as a function of the center of mass energy Ecu- The SMEC calculations have been done 
for the DDSMl residual interaction (^. For more details see the description in the text. 
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FIG. 9. The phase-shifts for the p+ O elastic scattering as a function of the proton energy Ep 
for the 1^, I"*", and ^ > | > | partial waves. The experimental data are from |^. SMEC 
results have been obtained for the DDSMl residual interaction and are plotted with the solid 
line. For I"*" partial wave, we show also the SMEC results without the resonant part (the dashed 
line). For more information see the description in the text. 
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FIG. 10. The elastic cross-section at a laboratory angle 9lab = 166° for the p + ^^O scattering 
as a function of proton energy Ep. The SMEC calculations (the solid line) have been performed 
using the DDSMl residual interaction. The results shown with the dashed line do not include the 
resonant contribution of 5/2^ state. Experimental cross-sections are from Ref. [40|. 
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